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We introduce a formula for the disagreement function which is used to 
control a recently proposed dynamics of the Ising spin system. This leads to 
four different phases of the Ising spin chain in a zero temperature. One of 
these phases is doubly degenerated (anti- and ferromagnetic states are equally 
probable). On the borders between the phases two types of transitions are 
observed: infinite degeneration and instability lines. The relaxation of the 
system depends strongly on the phase. 



1 Introduction 

The Ising spin system is one of the most frequently used models of statistical 
mechanics. Its simplicity (binary variables) makes it appealing to researchers 
from other branches of science including biology sociology and econ- 
omy 0, ^. In sociophysics models of opinion formation based on the social 
impact theory (reviewed in [||), the individual opinion is decribed by the 
Ising spin. This corresponds not only to typical "yes"-"no" questions, but 
also to important issues where the distribution of opinion seems to be bi- 
modal, peaked on extreme values. In general, in these models the influence 
flows inward from the border to the center, like in the majority rules, where 
the site in the middle takes the state of the majotiry of neighbouring sites. 
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In contrast, in our (USDF) model (reviewed in ^) an outward flow of 
influence is imposed. In the USDF model, an isolated person does not con- 
vince others; however, a group of people sharing the same opinion influences 
their neighbors. In spite of simple rules the model exhibited complicated 
dynamics in one p and more (references in |^) dimensions. In less than a 
year, this model has found several applications: e.g. it was used to explain 
the distribution of votes among candidates in Brazilian local election |^ and 
to model the price dynamics 0. 

In this paper we introduce the "disagreement function" which is used 
control the dynamics of the model. We show that for a one dimensional 
Ising spin chain in zero temperature this leads to four different phases: ferro- 
magnetic, antiferromagnetic, (2,2) antiphase and a doubly degenerated phase 
in which both ferromagnet and antiferromagnet are equally probable stable 
steady states of the system. Apart of structural differences between phases 
the difference in relaxation will be shown. The system in general will relax 
in two different ways depending on the phase. Moreover, a sharp change of 
the relaxation time on borders of the phases will be observed. 



2 The model 

Recently a simple model for opinion evolution in closed community was pro- 
posed P]. In this model the community is represented by a horizontal chain 
of Ising spins, which are either up or down. A pair of parallel neighbors 
forces its two neighbors to have the same orientation (in random sequential 
updating), while for an antiparallel pair, the left neighbor takes the orienta- 
tion ot the right part of the pair, and left neighbor follows the right part of 
the pair. Thus the model can be described by two simple dynamic rules: 

• Dv. Si.iit + 1) = S,it) and 5,+2(t + 1) = ^^(t) if S^it) * = 1 

. D2: 5,_i(t + 1) = 5,+i(t) and S,+2{t + 1) = S,{t) if S,{t) * = -1. 



In contrast to usual majority rules [jlO|, in this model the influence does not 



flow inward from the surrounding neighbors to the center site, but spreads 
outward from the center to the neighbors. The model thus describes the 
spread of opinions. The dynamic rules leads to two different stable steady 
states (ferromagnetic and antiferromagnetic) with equal probability. The 
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second dynamic rule (-D2) of the model has been already changed in two 
different ways. In case of antiparallel spins the neighboring spins can either 
flip with probability 1/2 |^ (-D2A) or remain unchanged [0 {D2b)- In both 
cases {D2A and D2B) the only final state is ferromagnet. It is worth to 
mention that the ferromagnetic state for both rules, D2A and D2B) is always 
reached (even in two dimensions) in contrast to the Ising spin system under 
Glauber dynamics |]lT], |12|. In the case of D2B besides of ferromagnetic stable 
steady states, the antiferromagnetic unstable steady state exsits. 

Since, we have up till now three different rules for the case of antiparallel 
spins, we propose a generalization of the previous models. The generalized 
model consits of two components, hence the name TC model: 

• The dynamics: choose a pair of spins S'j+i and Si^2 and change its next 
nearest neighbors Si and S'j+s 

• The rules: control the dynamics of the i-th and (z + 3)-th spins by the 
disagreement function. 

In the next sections we introduce the disagreement function and show that 
TC model includes as a special cases all earlier proposed models |, 0]. 
Moreover, TC model consits of more then those three subcases which we 
present on its phase diagram. Using Monte Carlo simulations we show how 
the system described by TC model relax. 

3 How to control dynamics? 

Let us assume for a while that we have the formula for a function that can 
control TC dynamics and denote it by E. We choose at random a pair of spins 
Si+i and 5*4+2 and we calculate = E{Si, Si-^-i, Si+2)- Next we calculate 
E^ = E{—Si, Si+i, Si+2) in the case of flipped i-th spin. If E^ < E^ then 
we will flip the z-th spin, if not the spin will remain unchanged. We do the 
same for the second neighbor of the choosen pair i.e. for the spin Si^^. 

Our dynamics looks now similar to the Glauber dynamics in zero tem- 
perature, where E plays the role of energy. However, there are three main 
differences between these two dynamics: 

• In Glauber dynamics we flip the i-th spin according to the interactions 
with (i — l)-th and {i + l)-th spins, here we look at {i + l)-th and 
(z -|- 2)-th spins. 
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In Glauber dynamics the flip is done even if the old energy is equal 
to the new one. It doesn't seem natural in zero temperature but it 
is needed to get the ground state (in two dimensions even this is not 
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• In our case E is called the disagreement function, since it is not the 
energy. 

Now we will look for the formula for E. We shall deal with the lattice model 
where each lattice site i is occupied by an ising spin Si = ±1. Usually, the 
spins are assumed to interact through pairwise coupling of the form —JijSiSj, 
where Jij are exchange integrals. Of course, the ordering of the spins is 
determined by the interactions. One of the best studied examples is the 
nearest neighbour (nn) Ising model with ferromagnetic coupling, i.e. Jij = 
J > for neighbour spins Si and Sj, while Jij = for more distant spins. 
Certainly, in a such model, the spins form the ferromagnetic state (all spins 
up or all spins down) in low temperature. For J < the anty ferromagnetic 
state is formed in low temperture. 

In TC model the i-th spin interacts with its two neighbors, and the ID 
hamiltonian can be written in the following form: 

H = —Ji^^ SiSi-^-i — J2^^ SiSi+2- (1) 

i i 

For Ji > and J2 < this is the well known ANNNI (axial next-nearest 
neighbour Ising) model introduced in [|13| and reviewed in It describes 



the Ising spin chain with ferromagnetic interaction Ji > between near- 
est neighbours (nn) and antiferromagnetic interactions between next nearest 
neighbours (nnn). Of course, in the one-dimensional case truly ordered states 
are stable only in zero temperature T = 0. If we introduce the competition 
ratio r = —J2/J1 we get in T = ferromagnetic state for r < 1/2 and (2,2) 
structure for r > 1/2. 

Now, we will use the nnn Ising hamiltonian ||l|| to construct the dissagre- 
ment function E. Chowdbury and Stauffer introduced similary a dissagre- 
ment function based on the simple nn Ising hamiltonian to the model of 
finantial market 0. We write E in the following form: 

E = —JiSiSi^i — J2SiSi^2- (2) 
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Each individual would like to minimize the corresponding disagreement func- 
tion. In the TC dynamics we choose a pair Si+i and S'j+2 and we change 
its neighbour Si (we also change S'j+3 spin calculating E = —JiSi+3Si+2 — 
J2Si+3Si^i, but for simplicity we further write only about i-th spin). For 
these three spins {Si, S'j+i, Si+2) we have four values of E: 

1. + + + or gives Ei = -(Ji + J2) 

2. - + + or H gives E2 = Ji + J2 

3. H h or - + - gives E^ = Ji - J2 

4. h or + + - gives E^ = J2 ~ Ji 

It is worth to notice that the possible transitions are only between states 
1 and 2 or between 3 and 4. Now we can derive from the TC model all 
previous models: 

• USDF model 0: 

if Si+i{t) * Si+2{t) = 1 then Siit + 1) = S,+i{t) i.e. Ei < E2 

if Si+iit) * Si+2it) = -1 then Siit + 1) = Si+2it) i.e. E3 < E,. 

Thus USDF model correspond to the TC model with — J2 < Ji < J2. 

• The model of financial market [^]: 

if * 5,+2(t) = 1 then S,{t + 1) = Si+iit) 

if Si+i{t) * Si+2{t) = -1 then Si{t + 1) = -Si{t) with probability 1/2. 
This corresponds to the TC model with Ei < E2 and E4 < E^ ^ 
—J2 < Ji and Ji > J2. 

• Other models reviewed in |]^: 

if S,+,{t) * S,+2{t) = 1 then S,{t + 1) = Si+,{t) i.e. E, < E2 
if Si+iit) * Si+2{t) = -1 then Si{t + 1) = Si{t) i.e. E^ = E^ 
These models correspond to the TC model with Ji = J2- 

There are of course more subcases of the TC model depending on interaction 
coefficients Ji and J2. On the Figure 1 all possible phases, depending on in- 
teraction coefficients, are presented. The North (doubly degenerated) phase 
corresponds to the original rule D2. The East (ferromagnetic) phase corre- 
sponds to rule D2A (the flip in case of antiparallel spins is made at random) . 
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The line between these two phases corresponds to rule D2B (the flip is possi- 
ble only in case of parallel spins). On this line the antiferromagnetic steady 
state still exists but it becomes unstable and we never reach it outside of this 
state. It is also interesting to see what happens on other border lines. The 
border between the ferromagnetic state and the (2,2) antiphase (see Fig. 2) 
is infinitely degenerated. Let us define (after [Q) a fc-band formed by k 
adjacent, identically oriented spins, terminated at the both ends by oppo- 
site oriented spins. With such a definition, the ferromagnetic structure is 
zero-band, antiferromagnetic is one-band and (2,2) antiphase is a two-band 
structure. On the line between ferromagnet and (2,2) antiphase any sequence 
of /c-band {k > 2) is equally probable (see Fig. 2). The line between (2,2) an- 
tiphase and antiferromagnet is also degenerated, and any sequence of k-band 
(with k=l,2) is the steady state (see Fig. 2). 

There is also another interesting feature which differs phases from each 
other - the time and the style in which the system relax. We will describe it 
in the next section. 

4 How does the system relax? 

What happens when we suddenly cool our system from a high temperature 
to zero temperature? As we mentioned previously the system will relax to 
one of the possible final states described by the phase diagram (Fig.l). But 
how does it relax? We studied this using Monte Carlo simulations. We found 
out that the relaxation process strongly depends on phase. The system can 
reach antiferromagnetic state in the West (antiferromagnetic) phase as well 
as in the North (degenerated) phase. However, it will relax to this state 
differently in each case. In the antiferromagnetic phase the system will be 
almost totally ordered after several Monte Carlo Steps (MCS). Then the 
system will oscillate oround the final state. These oscillations will decrease 
in time and finally the system will reach the steady state. In the degenerated 
phase the system will order very slowly. 

In Figure 3 the examples of relaxations in all four phases are presented. 
To show this relaxation we choose the opinion changes, since the model was 
proposed to investigate the opinion dynamics. We defined the opinion [^] as 
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a magnetization of the system: 

N 

m = J2Si, (3) 

1=1 

For such a choice the system will relax to \m\ = 1 (ferromagnet) or \m\ = 
(antiferromagnet or (2,2) antiphase). Of course, one could also choose the 
two point correlation function g =< SiSi^i > to see how the system relaxes. 
We have done it to recognize the final state ((7 = 1, —1 or for ferromagnet, 
antiferromagnet and (2,2) antiphase respectively). For Ji > —J2 (North 
and East part of the diagram in Fig.l) the ordering of the system is very 
slow. Sometimes the opinion can change dramaticaly in a short time (see 
Fig. 3). The long time trends are observed, which reminds very much of the 
real sociological procesess 0. For Ji < —J2 the system is almost ordered 
after several Monte Carlo steps, however, then it takes a long time to reach 
the real final steady state. The opinion is fluctuating around zero and these 
fluctuations are decreasing in time (see Fig.3). Although the way in which 
the system relaxes in the North and East phases is the same, the relaxation 
time in each of these phases is different. About a two times shorter (in 
average) time is needed to reach the final state in the degenerated phase. 
The relaxation time changes very sharply on the border between these two 
phases (Fig. 4). A similar effect is observed also on the border between the 
antiferromagnetic and degenerated phases. 

5 Summary 

We proposed the new generalized model of opinion formation. The disagre- 
ment function was introduced to control the simple dynamics of an Ising spin 
chain in zero temperature. This allowed to generalize the previous model of 
opinion dynamics. It was shown that the phase diagram for that system de- 
scribed by such a model consists of four different phases. The most interest- 
ing is the existence of the doubly degenerated phase in which the system can 
reach the antiferromagnetic steady state or the ferromagnetic steady state 
with the same probability. Moreover, it was shown that the system can relax 
in two different ways depending on the interaction coefficients. Surprisingly 
the system can reach the antiferromagnetic state in two different ways. In 
the antiferromagnetic phase the system will be almost ordered after several 
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Monte Carlo steps and then decreasing oscillations around the final state will 
lead the system into this state. In the degenerated phase, the system will 
behave "blindly" making a long "random" walk to the final state. It would 
be probably worth to look at the system described by such a model in higher 
dimensions and higher temperature. We also hope that the generalized TC 
model will find so many applications as its older brothers 0. 

I would like to thank the Fundation for Polish Science (FNP) for the 
finantial support. 
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Figure 1: The phase diagram of the TC model 
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(2,2) antiphase 




antiferro-antiferro degeneration (k=1 ,2) 



antiferro-ferro degeneration (k>1) 

Figure 2: Examples of Sdifferent steady states of the TC model are presented. 
Bright lines denote spins up and dark lines denote spins down. 
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(a) J^=0.8, J2=1 



(b)J^=1.2,J2=1 



time [10* MCS] 



time [10'* IVICS] 



E 




tlme[10'* IVICS] 



E 



(d) J^=-1,J2=0.5 



time [lO'' IVICS] 



Figure 3: Examples of the relaxation for 1000 spins system are presented. 
Two kinds of relaxations where observed depending on interaction coeffi- 
cients. For Ji > —J2 the system makes long "random" walk to the final 
state, while for Ji < — J2 the system makes decreasing oscillations around 
the final state. 
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J1 

Figure 4: Relaxation time for J2 = I. On this figure we present results for 
the system of 1000 spins averaged over 10000 samples. 
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